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ABSTRACT 

In this paper we consider two outstanding intertwined problems in modern 
high-energy astrophysics: (1) the vertical thermal structure of an optically thick 
accretion disk heated by the dissipation of magnetohydrodynamic (MHD) tur- 
bulence driven by the magneto-rotational instability (MRI), and (2) determining 
the fraction of the accretion power released in the corona above the disk. For sim- 
plicity, we consider a gas-pressure-dominated disk and assume a constant opacity. 
We argue that the local turbulent dissipation rate due to the disruption of MRI 
channel flows by secondary parasitic instabilities should be uniform across most 
of the disk, almost up to the disk photosphere. We then obtain a self-consistent 
analytical solution for the vertical thermal structure of the disk, governed by the 
balance between the heating by MRI turbulence and the cooling by radiative 
diffusion. Next, we argue that the coronal power fraction is determined by the 
competition between the Parker instability, viewed as a parasitic instability feed- 
ing off of MRI channel flows, and other parasitic instabilities. We show that the 
Parker instability inevitably becomes important near the disk surface, leading to 
a certain lower limit on the coronal power. While most of the analysis in this 
paper focuses on the case of a disk threaded by an externally imposed vertical 
magnetic field, we also discuss the zero-net-flux case, in which the magnetic field 
is produced by the MRI dynamo itself, and show that most of our arguments and 
conclusions should be valid in this case as well. 



Subject headings: accretion, accretion disks — magnetic fields — radiative transfer — 
galaxies: active — X-rays: binaries 
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Introduction 



1.1. Motivation 



Angular-momentum transport (AMT) in accretion disks is an important topic in mod- 
ern high-energy astrophysics, with applications to systems as diverse as Young Stellar Ob- 
jects (YSOs), accreting stellar-mass compact objects (white dwarfs, neutron stars, and 
black holes) in galactic binary systems (including X-ray binaries, XRBs), supermassive 
black holes in active galactic nuclei (AGNs), and even collapsar central engines of long 
gamma-ray bursts (GRBs). The leading candidate for explaining the observed relatively 
high AMT levels is the magn etohydrodynamic (MHD) turbulence excited by the magneto 



1991 



rotat i onal instabil ity (MRI) ( iVelikhovi Il959l ; I Chandr asekharl Il96ll ; iBalbus fe Hawley 
19981 ; lBalbusll2003l ). Over the past two decades, there have been many analytical and numer 
ical s tudies of MRI-turbulent disks [including three-dimensional (3D) MHD simulations, see 



g.. ^Brandenburg et al. 



1995 



1996 



Stone et al.lll996l: lArmitagJl998l: 



Miller fc Stone 



2000 



Turner 2004; ISano et al. 12004: Hirose et al. 



Davis et al. 2010 ; Lesur fc Longarettil 



2006 



Fromang fc Papaloizoul 120071; iFromang 



2010 



2011 



Beckwith et al. 



2009 



2011 



Simon et al. 



2009 



2011 



Shi et al 



2010 



Guan fc Gammie 



2011 



Bodo et al. 



2011 



Blaes et al.ll201lh ] re 



suiting in a significant progress in ou r understanding of how these systems work (see, e.g., 
Balbus fc Hawleylll998l ; IBalbus! 12003k for review). However, most of these studies have fo- 
cused, at least until recently, on the basic dynamical behavior of MRI-driven MHD turbulence 
in disks, with the primary goal of determining the rip component of the turbulent Reynolds 
and Maxwell stresses that drive accretion. Relatively little attention, in our view, has been 
paid to the questions of thermodynamics, in particular, to the vertical thermal structure of 
MRI-powered disks, governed by the balance between the heating due to turbulent dissi- 
pation (which is inextricably linked to the AMT) and the cooling by radiative processes. 
Correspondingly, most accretion disks simulations have focused on the MHD aspects of the 
problem, while treating the thermodynamics in a relatively simplified fashion, e.g., consider- 
ing an isothermal disk, or including optically thin radiative cooling. While these approaches 
have been very helpful in establishing the basic picture of AMT inside disks, they could not 
address the structure of the upper layers of disks, where the observable radiation is formed, 
nor the transition to the magnetically -dominat ed disk coron a. One notable excep t ion is the 
numer i cal work by 



5irose et all (120061 ) [see also Turner! ( 12004] ); iHirose et all (120091 ) : IShi et al. 



( 120101 ); iBlaes et al.l ( 1201 lh ] . who have made the first difficult strides towards a more real- 
istic treatment of MRI-powered radiatively cooled accretion disks. One, of course, has to 
appreciate the tremendous computational challenge of combining high-resolution 3D MHD 
simulations, necessary to get the basic magnetohydrodynamics of the problem right, with 
radiative transfer that spans both optically thick and optically thin regimes. 
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Another outstanding problem in today's high-energy astrophysics is understanding the 
formation of a strongly-magnetized hot tenuous corona above a turbulent accretion disk, in 
particular, understanding what governs the fraction of the accretion power that is released 
in the corona above a relatively cold and dense disk. The observational motivation for our 
interest in accretion disk coronae (ADCe) stems from the fact that they lie at the base of 
disk-driven outflows and thus link disks to their winds and jets and also from their role as 
emitters of high-energy (e.g., X-ray) radiation in many different types of accreting systems. 
Examples of astrophysical systems were the presence of an ADC has been inferred include : 
galactic black-hole (BH) XRBs such as Cyg X-l fe.g.. iBisnovatvi-Kogan fc Blinnikoyl 11976 



Liang fc Pricdll977l ; iGierlinski et al. 



mass ive BHs in AGNs (IHaardt fc Maraschi 



1999 



C 



lurazov et al 



1991 



1993 



2001 



Di Matteo et a,. 



1982 



Kawaguchi et al 



2001 



1999h: super - 



Kriss et al. 



1999 ) : and accreting white-dwarf binary systems (Cataclysmic Var iables, CVs) (e.g.. lWhite fc Holt 



Church et al. 



19981 : iNavlor et al.lll988l : iRamsav et all 120011 ): and YSOs, e.g., T Tauri 



stars, (IKwanl 119971 ) (the disk coronal power in these systems is likely to manifest not in 
X-rays but at lower energies, e.g., UV and optical). 

In studying the formation of ADCe, one would like to understand how much magnetic 
flux and associated magnetic energy emerge from the disk into the overlying corona and what 
physical processes control the rate and the form of the buoyant magnetic flux emergence. 
An issue of particular observational importance is the significant variation in the observed 
relative levels of (coronal) X-ray activity among different spectral states in galactic black- 
hole X-ray binaries and also among the different types of black-hole accreting systems. In 
particular, the coronal fraction in AGNs can often be as high as a few tens of percent (e.g., 
Svensson &: Zdziarskilll994l : IWang et al.ll2004j ). significantly higher than in galactic black-hole 
X-ray binaries in the high-soft state. Understanding the reasons for this dichotomy probably 
requires a detailed quantitative picture of MHD turbulence in a stratified radiation-pressure- 
dominated disk (which is, after all, the ultimate energy source for coro nal activity), with an 
emph asis on the production and buoyant rise of magnetic structures (IBlackman &: Pessah 
2009h . 



A major theoretic al development in our understanding of ADC formation was the model 
of iGaleev et al.l ( 119791 ). who described the key general magnetic processes leading to the 
em ergence of the corona. Our conceptual theoretical understanding was further advanced 
by iTout fc Pringld (119921 ) who outlined a dynamo cycle based on the interplay of the MRI 
and Parker instabilities. 



-4- 



1.2. Comments on the Physics of MRI Turbulence 



In local shearing-box numerical studies of MRI turbulence, it is customary to distinguish 
two cases: the case with a finite net vertical magnetic flux, which i s conserved in the c ourse 
of the simulation, and the case with a zero net vertical flux (e.g., lHawley et al.lll995[ ). Al- 



though both cases have been studied extensively in the past, in recent ye ars it seems that the 



focus has shifted towards the zero net flux case (except for the work by iLesur fc Longaretti 



201ll ). We note, however, that the non-zero vertical flux case is still of considerable interest, 
especially in situations where one is interested in the interaction between the turbulent ac- 
cretion disk and the large-scale magnetosphere. In particular, one may be interested in the 



bulent disk ( 


van Balleeooiier 


i 1989 


Lubow et al. 


1994 


Livio et al. 


1999; 


Spruit & Uzdensky 


2005; 


Rothstein & Lovelace 


2008 


: ] 


3eckwith et al 


J 2009 


). One of the most notable examples 



of such a situation is the magnetic interaction between the disk and the dipole-like magnetic 
field of the ce ntral star, e.g., in the context of Young Stellar Objects and Neutron Star X-ray 
binaries (e.g., lUzdenskyl 120041 ) . Another example is when one is interested in a large-scale 
unipolar magnetic field that is believed to be important for launching disk outflows (winds 
and jets) in YSOs, galactic BH XRBs, AGNs, and collapsar central engines of GRBs. 

In accretion disks threaded by a relatively weak (below equipartition with the gas pres- 
sure) vertical magnetic field with a non-zero net flux, gravitational and rotational energy of 



the accreting matter is first transformed by the MR 



the so-called chan nel flows (lHawley fc Balbus 



1992 



into the magnetic and kinetic energy of 



Sano fc Inutsuka 2001; Sano et al. 2004 



Bodo et al.ll2008[ ). As was shown by lGoodman fc Xul (119941 . hereafter GX94) in the shearing 
box formulation, these channel flows are exact nonlinear solutions of the MHD equations; 
they are basically MRI linear eigenmodes that continue to grow exponentially even when 
they become very large (IGoodman fc Xulll994l ). However, as was also shown by GX94, in 
reality, channel flows cannot grow indefinitely, since they themselves become unstable to 
various secondary parasitic instabilities — i.e., instabilities that feed off of the velocity and 
magnetic field structures produced by the primary MRI mode, which in this context plays 
the role of a slowly evolving equilibrium. The development of the parasitic modes leads to 
the disruption of MRI channel flows and a fully-developed turbulence that dissipates energy 
locally in the disk. GX94 themselves considered only ideal-MHD parasitic instabilties in the 
non-stratified case and found two families, one of which is related to the Kelvin- Helmholtz 
instability. In principle, however, other parasi tic modes are possible, such as non-ideal , visco - 
resistive modes including the tearing mode ( iPessah &; Goodman! |2009| ; lLatter et al.ll2009l ). 
It is also important to note that magnetic fields generated by the MRI channel flows are 
mostly toroidal and are themselves subject to a toroidal-field version of the MRI, which thus 
should be viewed here as a secondary parasitic instability feeding off the primary MRI mode. 
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Another limitation of the above-mentioned analytical studies of parasitic modes is that 
they restricted themselves to considering an unstratified shearing box; this can be thought 
of as a representation of a small region inside the disk, small compared with the pressure 
scale-height, H ~ c s /Q. Thus, these analyses were effectively local not only in cylindrical 
radius but also in the height z above the disk midplane. Since the characteristic scale of 
the fastest-growing MRI mode, / mr j, is proportional to the strength of the vertical magnetic 
field, 

U = ^~vk/n, (l) 

Z7T 

where Va is the Alfven velocity corresponding to the local density and the vertical magnetic 
field B , this approach is justified only if the vertical magnetic field is sufficiently weak, i.e., 
Va "C c s . The generalizatio n of the channel mo de analysis to the case of stratified disks has 



been developed recently by lLatter et al.l (120101 ) 



We would like to remark that recently reported numerical evidence from large shearing 
box simulations suggesting t hat channel flows are not im portant in the overall energetics and 



AMT in accretion disks (e.g.. lLongaretti fc Lesurll2010l ) can probably be attributed to a very 



restricted definition of channel flows used in these studies, in particular, to restricting it only 
to axisymmetric modes in the turbulent Fourier spectrum. Of course, if the computational 
box is sufficiently large in the toroidal direction, then the contribution from axisymmetric 
(i.e., spanning the entire toroidal extent of the box) modes should indeed be small. This is 
because any coherent flow and magnetic field structures, such as channel flows, have only a 
limited lifetime before they are disrupted by the parasitic modes as discussed above and thus 
cannot be correlated over distances much larger than about the MRI wavelength. That is, 
one should not realistically expect channel-mode structures extending (in any direction) to 
sizes much larger than A mr i. Since here we are interested in a situation where A mr i <C H, then 
the contribution of any axisymmetric (or, for that matter, any large-scale, / ~ H) structures 
in the overall energy dissipation and AMT should indeed be inevitably small. We believe, 
however, that this does not disqualify any localized channel-mode structures. Locally, on 
scales of order / mri <C H, we can still expect that the picture presented in GX94 basically 
holds. 



1.3. Objectives of This Paper 

The present paper has two main goals. First, we are interested in the vertical structure 
of a gravitationally stratified accretion disk. For definiteness, we will restrict the present 
study to the case of a gas-pressure dominated disk, even though it is not directly applicable 
to inner parts of black-hole accretion disks in the high soft state. Developing a generalization 
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of our model to the case of a radiation-pressure-dominated disk is left for a future study. 

As ment ioned above, most of the previous numerica l studies of MRI turbulence (with the 
exception of iHirose et al.l 12006k 120091 ; iBlaes et al.l 1201 ll ) were done using either the isother- 
mal approximation or a simple prescription for optically thin radiative cooling. Standard 
geometrically thin accretion disks, however, are optically thick, and so one inevitably has 
to face the optically-thick radiative transfer problem in order to deduce their structure 
( IShakura fc Sunyaevi 119731 1. Thus, the first main goal of the present study is to construct a 
self-consistent model of the vertical structure of an optically thick accretion disk heated by 
the local dissipation of MHD turbulence and cooled by radiative diffusion (we shall ignore 
any external irradiation). To attack this problem, for definiteness we will adopt the point of 
view that the path to turbulence onset lies mostly in the nonlinear disruption of the MRI 
channel modes by parasitic instabilities, as described above. This assumption will allow us to 
obtain a concrete vertical profile (namely, flat, see § 13. ip of the magnetic energy dissipation. 
We note, however, that this assumption is probably not critical and a similar picture may be 
developed without an explicit reliance on the concept of parasitic modes. In any case, having 
a specific physically motivated prediction for turbulent dissipation, supplemented with some 
extra assumptions about the radiative transportproperties of the disk (see §E]), will allow 
us to construct a full t heory of the steady-statqj vertical disk structure, along the lines of 



the classical analysis of IShakura fc Sunyaevi ( 119731 ) 



Our second main objective is to find a way to estimate the fraction of the accretion 
energy released by the MRI that is not dissipated locally in the disk itself, but instead is 
transported vertically by buoyantly rising magnetic flux tubes and dissipated in the hot 
overlying corona. In our view, the coronal fraction of the released power is, to a large de- 
gree, deter mined by the competion between the GX94 parasitic instabilities and the Parker 
instability (lParken\l96a) , viewed here as another t ype of parasitic instability feeding off the 
horizontal magnetic field of the p rimary MRI mode ( jFoglizzo &; Taggerlll995l ; iTout fc Pringle 
19921 ; iBlackman fc Pessahl 120091 ). As a first step towards investigating this issue quantita- 
tively, we will compare the growth rates of the GX94 and Parker parasitic instabilities and 
will investigate how their ratio varies with height. As we will show in this paper, the Parker 
instability is slower than other parasitic instabilities in most of the bulk of the disk, but 
starts to become competitive in the disk's upper layers, leading to a certain lower limit on 
the coronal power. 

In most of the paper, we will focus on the case where the main magnetic field respon- 
sible for the MRI is an externally imposed large-scale vertical magnetic field B , uniform 



1 We are here interested in intermediate timcscales At such that 1 , Ai ra d diff < < At a 
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in z. This is applicable to situations where such a field is relatively strong, namely, stronger 
than the magnetic field that would be produced in its absence by turbulent dynamo ac- 
tion due to the MRI turbulence itself. We note that, in general, the problem of large-scale 
dynamo in stratified MRI turbulen t disks is still poorly understood and represents a key fron- 



tier i n the accretion disk research ([Brandenburg et al. 



1996 



2009; Simon et al. 



Balbus &: Hawlev I1998L iBlackman &: Tanl 12004 ; 



1995 



Stone et al. 



Vishniac 



2009 



1996 



Hawlev et al. 



Blackman &; Pessah 



20121). Nevertheless, as most of t he recent zero-net-flux numerical sim- 



ulations indicate ( iShi et al.l |2010| ; iDavis et al.l |2010| ; iGuan &; Gammid 120111 ). the MRI dy 



namo, by itself, leads to a saturation large-scale magnetic field only at a relatively low level, 
-Bdyn *C -B eq = (^irPo) 1 ^ 2 , where P is the midplane gas pressure. Moreover, even the small- 
scale turbulent magnetic field responsible for the A MT in zero-mean flux simulations is also 
relatively small, of order I? t 2 urb /87r ~ 0.01 Pq (e.g., IBlackman et al.ll2008l ). Thus, if the ex- 



ternal mean vertical field B is larger than this characteristic dynamo field but still weak 



compared to B eq , i.e. 



2 

turb 



« B 2 « B 2 then one may expect the dynamics to be 



determined mostly by the super- impo sed vertical field Bn . This point o f view is also sup- 
ported by recent numerical studies by iLesur fc Longarettil (120111 ) and by iBodo et al.l (120111 ) 
and we will adopt it as our starting point in the main part of the paper where we consider 
the non-zero net flux case. 

We note, however, that the opposite case — the case where the externally imposed 
magnetic field is small or absent and where the MRI is driven by the self-generated dynamo 
magnetic field, - - is also, of course, of considerable interest in astrophysics and we will 
devote § |7] to discussing it. In particular, we will argue that most of our results for the 
vertical disk structure obtained in § [3] for the finite- B case can also be applied to the zero- 
net-flux case , with Bp replaced by -Bhv ti- Furthermore, we will recover the basic scalings 
obtained by IShakura fc Sunyaevl (119731 ) for this case and, in addition, will calculate the 
vertical temperature and de nsity profiles t hat se em to be in good agreement with the results 
of numerical simulations by iHirose et al.l (120061 ) . However, we will not attempt to estimate 
the coronal power fraction for the zero-net-flux case, leaving this task for a future study. 

This paper is organized as follows. In § [2] we describe the basic idea and our overall 
approach to the problem. In § [2] we present the calculation of the internal vertical thermal 
structure of the disk threaded by a nonzero net vertical magnetic flux. Next, in § H] we 
consider the upper layers of the disk and investigate how close to the disk surface various 
assumptions of our model break down. Then, in § [5] we estimate the growth rate of the 
Parker instability and compare it with that of the GX parasitic instabilities, as a function 
of height; this allows to estimate the coronal fraction of the accretion power as a function 
of the system's parameters. Then, in § [6] we discuss the implications of our model for the 
accretion rate and for the longer-term evolution of the accretion disk. The zero-net-flux case 
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is discussed in § [7J Finally, we present our conclusions in § [HJ 



2. The Overall Physical Picture 

The overall physical picture we have in mind can be described as follows. 

For MRI to be active in the first place, the vertical magnetic field must be relatively 
weak, with the corresponding magnetic pressure less than the gas pressure P in most of 
the disk, Bq/8tt < P. [As mentioned above, for simplicity we ig nore radiation pressure 



in the present work.] Furthermore, followin g the previous studies ( IGoodman fc Xul 11994 



Pessah &: Goodmanll2009l ; lLatter et al.l 120091 ). for simplicity we shall consider the case when 
this field pressure is not only weak, but very weak compared with the gas pressure: 5q/8/T <C 
P. Having the small parameter = Bq/8ttP <C 1 at our disposal will allow us to make 
several important simplifications. 



First, in the spirit of GX94 and lPessah fc Goodman! ( 120091 ). we assume that the growth 



of MRI is checked by the development of the parasitic instabilities. This assumption allows us 
to get on a direct path towards evaluating the MRI turbulent energy dissipation ratej^l Once 
the parasitic instabilities take over, fully-developed MHD turbulence sets in and destroys the 
horizontal magnetic field of the primary MRI mode. As a result, the magnetic and kinetic 
energy of the primary MRI mode are dissipated by turbulent cascade on a timescale of order 
the dynamical time, i.e., Q" 1 . The typical maximum amplitude to which the MRI channel 
flows are able to grow can be represented by the characteristic horizontal magnetic field 
-Bhor,sat = b sat B . It can be estimated from the condition that the growth rate of the fastest- 
growing parasitic mode is comparable to the growth rate of the fastest growing (and hence 
the most relevant) primary MRI mode, which is essentially a numerical constant times Q, 
independent of (5. The maximum growth rate of the GX94 parasitic instabilities is naturally 
proportional to the amplitude of the primary mode, i.e., 'Jgx ~ bQ. Therefore, the two 
growth rates become equal at a certain finite critical value b = 6 sat , generally of order one 
(even though the GX94 theory was formally developed under the assumption that 6 > 1). 
Importantly, within the framework of the GX94 model (unstratified disk with /3 > 1, ideal- 
MHD incompressible motions, etc.), there are really no additional parameters on which the 
dimensionless saturation amplitude 6 sat could depend. It is then natural to take it to be just 



2 We note, however, that our model should not, in our view, rely critically on our assumptions about the 
role of parasitic modes in setting the saturation level of MRI turbulence. We believe that essentially similar 
results can be obtained invoking some other physical picture of the saturation mechanism, e.g., turbulent 
diffusivity (J. Goodman, private communication). 
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a constant number of order unity. In particular, it should be the same for all heights in the 
disk, as long as we are not close enough to the disk surface, where the assumptions of the 
model break down (see § H]). 

Another reason why the assumption Bq/8tt <C P is important is that it justifies the 
locality (in z) of our picture of MRI turbulence. Indeed, the communication speed in the 
vertical direction is only Va, and the lifetime of channel modes limited by the disruption 
by the parasitic instabilities is of order Q~ 1 \ log 6 sat I - Any causal connection established 
during this time extends only over a vertical distance of order / mri | log £» sat | ~ Z mri , which is 
smaller than the gas pressure scale height H by a factor /3 1 / 2 . This means that, as long as 
our assumption 3> 1 holds, channel flows that develop at substantially different heights 
interact with each other only weakly and hence can be considered separately. In other words, 
MRI at different heights in the disk develops independently of what happens at other heights. 
Therefore, since we are interested in the vertical structure of the disk including the vertical 
profile of the MRI turbulence and its energy dissipation rate, we can just regard properties of 
MRI channel flows and parasitic instabilities as being local in z. One thus can conclude that, 
under these circumstances, the exersize of calculating global MRI eigen-modes spanning the 
entire thickness of a stratified disk may be of purely academic interest. 

The dissipated turbulent energy provides the main heating source of the gas in the disk. 
Because MRI turbulence is distributed over the disk thickness, the effective heating source is 
also distributed, Q — Q(z). In a thermal steady state, this distributed heating is balanced by 
the sum of turbulent thermal conduction losses (which, as we shall argue below, effectively 
just modify Q by a constant factor of order unity) and radiat ive losses; we assume tha t 



the latter take place via optically thick radiative diffusion (e.g., IShakura fc Sunyaevlll973l ). 
Solving the corresponding vertical radiative transfer problem yields a full self-consistent 
vertical structure of the disk (see 



Having this solution at hand will help us obtain a lower estimate for the coronal power 
fraction (see §EJ). As mentioned in § 11.31 our approach to this problem is to regard the Parker 
instability, which leads to the buoyant rise of the magnetic flux tubes into the corona, as a 
parasitic instability that competes for power with the other (GX94) parasitic instabilities. 
As we shall show in this paper, deep inside the disk, at heights z < H ~ c s /Q, the Parker 
instability growth rate is smaller than the characteristic GX94 growth rate, roughly by a 
factor of /3 1 / 2 > 1. This means that in the bulk of the disk, the MRI channel flows are 
destroyed by the GX94 parasitic instabilities well before the Parker instability can develop. 
Correspondingly, most of the energy density associated with the MRI channel flows (~ 
£> 2 or /87r ~ b 2 E>l) goes into feeding the resulting local MHD turbulence and is dissipated 
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locally in the diskjfl However, in the upper layers of the disk, as one approaches the disk's 
photosphere, the Parker instability starts to compete effectively with the other parasitic 
instabilities and hence one can expect a significant fraction of the MRI energy in this region 
to be transported up into the corona. 

Whether the assumptions on which the above arguments are built are actually valid, 
will need to be investigated by carefully designed numerical simulations of MRI turbulence in 
a stratified accretion disk, including full gas thermodynamics with radiative transfer. Since 
the present paper provides specific predictions for the disk structure and the coronal power, 
we hope that it will motivate such studies in the future. 



3. Vertical Structure of a Disk Threaded by a Net Vertical Magnetic Field 

3.1. Vertical Profile of Energy Dissipation 

The first important conclusion we can derive from the arguments in the preceding sec- 
tion is that, as long as the parasitic instabilities dominate, the volumetric dissipation (i.e., 
heating) rate should be independent of height z within the disk. 

Indeed, the vertical magnetic field is independent of z, but the density decreases with 
height because of stratification. This means that Va increases with z, and so does l mri . 
Then, both the typical vertical and horizontal extent of MRI channel modes, of order Z mr ; 
and Z m ri&sat, respectively, increase with height; however, as long as they remain <C H, the 
typical life time of the channel modes, and the local turn-over timescale of the MRI turbulence 
remain constant, of order The typical horizontal magnetic field component also remains 
independent of height, of order B hoT ~ B b sat . The overall local dissipation rate per unit 
volume can then be estimated as 

q ~ 7mri 3« tog-^/aStetf))] ~ n bl t . (2) 

Here the logarithmic factor represents the number of e-foldings needed to grow from some 
initial perturbation 5Bh or (0) to the saturation amplitude -B sat . In a realistic situation, of 
course, one does not expect the magnetic field to return exactly to a pure vertical-field state; 
therefore, in the following analysis, we shall assume that the typical initial perturbation 
amplitude is of the order of Bq itself, and thus will ignore this logarithmic factor. 



3 We note, however, that a small fraction of this energy may still escape magnetically into the 
corona if there is an effic ie nt inverse-ca s cade dynamo producing large enough magnetic structures, e.g., 
Blackmail fc Pessahl (|2009h : ISimon et all (|2012h . 



- 11 - 



Thus we now see that, because the background vertical field threading the disk, B z = Bq, 
is independent of height z in the case under consideration, the above volumetric MRI dissi- 
pation rate is also independent of height, basically unaffected by the pressure stratification! 
Importantly, it does not have to follow the gas density or the gas pressure. That is, the 
volumetric dissipation rate does not follow mass, as it is sometimes assumed, but is just 
constant, uniform in z, at least across the main part of the disk, where the magnetic field is 
still dynami cally weak. This expectation is in fact c onsistent with the results of numerical 
simulations ( iMiller &: Stondl2000t iHirose et al.l 120061 ) for the zero- net-flux case, which show 
turbulent magnetic dissipation that is roughly flat over a couple gas scale-heights (or even 
slightly peaked at about z = 2H). 



It is important to note that, since the ultimate source of the heating is the dissipation 
of the accretion energy, the fact that Q ~ Bq implies that the accretion torque is also 
proportional to the square of the mean vertical field. This is in clear contradiction with 
the results of a number of numerical simulations rep orting a linear scali ng of the ac c retion 
stress with B (IHawley et al.lll995l ; |Pessah et al.ll2007l ) [note, however, that lSano et al.l (12004 ) 



report a B^ 2 scaling]. We believe that this discrepancy may be attributed to the fact that 
most of these numerical studies were not in the asymptotic regime of interest here, i.e. they 
lacked the required separation of scales between the disk scale height (or the vertical box 
size L), the prevailing MRI scale Z mr ;, and the dissipative scale (or the grid scale). In other 
words these studies either were not sufficiently resolved, in the sense that their MRI scale 
corresponding to the net vertical field B were not much larger than the resolution scale, 
or their MRI scale was not much lower than the scale height (for stratified sims) or the 
box size (for unstratified sims), which is equivalent to saying that Bq was not sufficiently 
small compared with the pressure-equipartition field (and hence the corresponding MRI 
waveleng th was not much smaller than the H). In fact, a recent very careful numerical 
study by iBodo et al.l (120111 ) demonstrates that indeed the accretion torque indeed scales as 



£>g in the asymptotic regime 5 <C l m r\ •C L, and hence s upports the point of view advanced 



in the present paper [see also lLongaretti fc Lesurl (120101 )] 



Another important point here is that the specific (i.e., per unit mass) AMT due to 
MRI turbulence and hence the resulting effective accretion inflow velocity are not uniform 
in height: at larger heights, the accreti on inward drift v elocit y is higher. This is similar to 
the results of numerical simulations of iBeckwith et al.l (120091 ) where this property of MRI 
turbulence in a stratified disk provided a mechanism for efficient inward transport of the 
vertical magnetic flux (see also iRothstein fc Lovelace] 120081 ). 
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3.2. Vertical Structure of an Optically Thick Disk 

Once the vertical profile of turbulent dissipation is est ablished, we can det e rmine the 



vertical structure of the disk, similar to the calculation by IShakura fc Sunyaevl (119731 ) for 
accretion disks and to standard stellar structure calculations. Here we are interested in a 
steady-state profiles of the gas temperature and density established on timescales longer than 
the characeteristic orbital time and the cooling time, but shorter than the overall accretion 
time (characteristic radial transport timescale). 

We also neglect external irradiation of the disk, e.g., by the radiation coming from the 
inner part of the disk or the central star. 

It is widely recognized that in discussing MRI in the presence of a mean vertical field, 
the ratio of the magnetic pressure to the gas pressure, f3~ 1 = Bq/StiP, is an important 
parameter. However, it is important to realize that the gas pressure profile in the disk is not 
just some arbitrary prescribed function, but needs to be determined self-consistently using 
the physical laws, namely the energy transport and the vertical hydrostatic balance. This 
point was made by Hirose et al. (2009), and here we adopt this point of view. That is, we 
cannot prescribe the temperature and density profiles of the disk; they are to be determined 
as part of the overall problem. Instead, we can prescribe as fixed only the those quantities 
that evolve only relatively slowly, namely, due to the radial transport associated with the 
accretion process itself; these quantites are then conserved on the timescales of interest to 
us here. In the present problem there are basically two such quantities: the disk surface 
density of mass, S, and that of the vertical magnetic flux, i.e., the vertical magnetic field B . 
In addition, we can prescribe the local disk rotation rate Q, and the parameters describing 
the radiative transfer, e.g., the scattering opacity n, which for simplicity we assume to be 
constant, see below. Thus, S, Bq, Q, and k are the only input parameters determining the 
disk vertical structure in our model. 

It will also be convenient to define a important dimensionless parameter — the disk's 
optical depth measured from infinity to the disk midplane: 

r tot = — , (3) 

We assume the disk to be optically thick, r tot 3> 1. 

Another important dimensionless parameter in our problem is the midplane plasma-/? 
parameter, defined with respect to the superimposed vertical field B : 



(4) 
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where Pq = P(z = 0). This parameter characterizes the dynamical importance of the vertical 
magnetic field, relative to the gas pressure. In our model we assume that /3 3> 1. Note 
that since we haven't yet computed the midplane plasma pressure, we cannot, at this stage, 
express (3q in terms of our principal input parameters; this will be done at the end of this 
section. 

Now let us solve for the vertical structure based on the hydrostatic bal ance and the radia- 



tive en ergy transport. This calculation is similar to the classical analysis by lShakura fc Sunyaev 



( 119731 1 and is also similar to the traditional analyses of the radiative stellar structure except 



that it is performed in plane, rather than spherical , geom etry. One important difference be- 



tween our analysis and that of IShakura fc Sunyaevl ( 119731 ) is that here we are able, within the 



framework of our model, to obtain explicit analytical expressions for the dependences of the 
plasma density p(z) and temperature T(z) on the height z within a gas-pressure-dominated 
disk. 

We shall start with the condition of hydrostatic pressure balance. Neglecting magnetic 
and radiation pressure support compared to the thermal gas pressure, and assuming that 
gravity is from the central point mass, we have (for a fully ionized hydrogen plasma): 

dP 2 d(pk B T) 2 

= 9zp = -M zp. (5) 



dz nip dz 



That is, 



4 d(pk B T) 
m p d(z 2 ) 



tt 2 p. (6) 



Next, we need to supplement this equation by the vertical heat balance equation that 
reads, in a steady state, 

Q — (-^rad + Fturb) ; (?) 

where _F ra d is the radiative energy flux and -F tur b is the effective vertical heat flux due to the 
MRI turbulence itself. The latter can be estimated roughly as a diffusive flux F tnr ^(z) = 
— D turh d(3nk B T)/dz, with an effective MRI-turbulent diffusion coefficient -D tU rb ~ 
Then, ignoring factors of order unity, we can write: F tuI ^(z) ~ — l mr iVA dP/dz ~ — Q^ 1 V\ dP/dz. 
Substituting dP/dz from equation (jSl), we get 

F Uah {z)~nv2z P = n?±z, (8) 

and hence, since both Q and Bq are constant in z, 

dF tmh Bq 

— - — ~ U — = const . 9 
dz Att 
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Thus, we see that the effective cooling rate due to the vertical heat transport by the MRI- 
driven turbulence is basically the same as the MRI turbulent heating rate Q (see eq. [2]), 
apart from a constant numerical coefficient of order unity: dF tnT \,/dz = ftarbQ] m particular, 
importantly, it is constant in z. This allows us to combine both of the effects of the MRI 
turbulence — the heating by turbulent dissipation and cooling by vertical turbulent transport 
— into one single term, the reduced heating rate 

dF 

Q' = Q 7— = riQ = const , (10) 

dz 

where we took into account our result that Q'(z) is constant (see § 13.11) and where we 
introduced a constant dimensionless factor rj = 1 — / tur b < 1 to account for the reduction, 
due to the turbulent heat transport, of the required radiative losses. 

Correspondingly, the energy balance equation now becomes 

dF iad /dz = Q' = const F rad (z) = Q' z , (11) 

In the radiation diffusion approximation, the vertical radiative energy flux is 

ca dT 4 ca dT A ca dT 4 

-Trad = — —j— = - — —7— \h = ~ ~ J" , (l^J 

6 dr 6 dz 3 pndz 

where a = 4<jsb/c = n 2 k\ 3 /l'bh z c 3 ' = 7.566 x 10~ 15 erg cm~ 3 K~ 4 is the radiation constant, 
a sb is the Stefan-Boltzmann constant, and A p h = —dz/dr = (p^)" 1 is the photon mean free 
path. This allows us to write our energy balance equation as 

2ca dT 4 



3 d{z 2 ) 

From the two equations (1131) and ([6]) we get: 



P kQ' . (13) 



^dT 4 = - KQ'pd{z 2 ) = ^ — d{pk B T) . (14) 

In general, the opacity should be a function of the local plasma parameters such as 
the temperature, but for simplicity and definiteness in this study we shall assume that it 
is constant, n = const (as it would be for the case of electron scattering, for example). 
Investigation of more realistic and complicated opacities, including the transition to the 
free-free opacity ~ pT~ 7 / 2 in the upper, colder layers of the disk, is left for future work. 



We can now immediately integrate equation ( 114]) to obtain an algebraic relationship 
between density and temperature: 

'M = S r3(2) ^ r3 ' (15) 
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where we defined 

= caQ 2 m p _ 2a SB Q 2 m p 

~ 6KQ'k B ~ 3KQ'k B ' [ ' 

and where we neglected the integration constant by assuming that T(z) and p(z) effectively 
become small together at the disk's photosphere. [We assume that the disk is not subject to 
any substantial external irradiation (such as that coming from the central star for example), 
which would invalidate the above boundary condition]. 

At this point we can check that a disk with this structure is always convectively stable 
for a gas adiabatic index equal to 7 ad = 5/3. Indeed, because p(z) ~ T 3 (z), it follows that 

V^=«) = i<V„ = I^4 (17) 

VCtlogP/rad 4 7 ad 5 

This means that the disk is stable to thermal convection and hence our assumptions regarding 
the nature of vertical heat transfer are justified. 

Substituting equation (ITBl) into (jHJ), we obtain: 

dT* _ m p tf 

I*" " ~ 1*7 dz ' (18) 

and hence 

k B T(z) = k B T(0) - ™^ z 2 = k B T(0) (l - 4) , (19) 

-LD V Zq / 

where T(0) is the midplane temperature and 

z = 4c sfi /Q (20) 

and c 2 Q = k B T(0)/m p . 

Next, using equation (I15|) . we obtain the density profile: 

p(z) = p(0) (l - = ^T 3 (0) (l - . (21) 
\ z / \ z Q / 

That is, in this model the disk has a sharp surface (similar to stars) at a finite height z — Zq 
above midplane. In the following, it will sometimes be convenient to normalize z by zq, i.e., 
to use the dimensionless height variable 

C = z/z . (22) 

The midplane values of temperature and density are then determined from the condition 
that the surface density has to be equal to a prescribed value: 



2 / p{z)dz = £. (23) 
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Thus we get 



2z oP (0) (l-CrdC = C lZo p(0) 



(24) 



where C\ = 2 /^(l — C 2 ) 3( ^C — 32/35. Of course, our model is too crude to consider any 
factors of order 1 meaningful, but we'll keep them here anyway. Thus, the total optical depth 
from infinity to the midplane of the disk and the central density are related via 



(25) 



Remembering the definition ( 1201) of zq, we then get the central temperature in terms of 
the primary input parameters: 



T = T(0) 



128 ft V 1 p! a SB 



2/7 



2/7 



105 B 2 bl tV 



128 8 



7T 



(k B /m p 



,V2 



r tot 
0"SB 



2/7 



(26) 



where we used Q' = r\ ft £> 2 & 2 at /87r as the fiducial reduced characteristic MRI heating rate, to 
get the last expression. It is interesting to note that the above expression for the midplane 
disk temperature involves only the vertical magnetic field pressure and the total disk optical 
depth, plus some universal physical constants. In particular, one can see that T(0) does 
not depend on the orbital frequency ft. This can be understood qualitatively as follows. 
The total dissipation rate per unit disk area is of order z Q f ~ (c Si0 /ft) (77ft B^b^/Sn) ~ 
T 1 ' 2 ^) 7] (i?g6 2 at /87r) - - independent of ft. On the other hand in this model this dissi- 
pated power is emitted from the disk's photosphere as thermal emission: F ~ <7sBTp h ~ 
cr SB^ 4 (0)r t ~ t 1 . Comparing these two expressions, we immediately find the above scaling of 
T(0) with Bl and r tot . 



Numerical estimate yields: 

T(0) 



220Kr t 2 o { 7 [7 7 S ^ at /87r 



12/7 



(27) 

Thus, in the case of galactic BH XRBs, for example, with a typical B ~ 10 7 G, we get 
T(0) ~ 10 6 K T^t, which is not unreasonable. However, in reality we, of course, do not 
expect our present model to apply in the inner parts of a black hole accretion disk because 
of the predominance of radiation pressure there. 

With the above result for the midplane disk temperature at hand, we can estimate other 
important disk quantities, e.g., the speed of sound 



c s,0 



j k B T(0) rl05 Q' 



128 ft 



k B 



Ttot 
0"SB 



1/7 



{21 
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and the disk half-thickness: 

z = Ac sfi /vt = 2ir 8/7 



105 



fcjA 4 QVtot j 1 / 7 
mj <T S B - 



(29) 



Next, the central density becomes 
po = p(0) = AT 



_ (35r tot ) 6/7 



The gas pressure profile is given by 



P(z) 



2pk B T 2Ak B 



T 4 (z) 



32 k 



2Ak 



3Q' \k B - 



1/7 



(30) 



Lp trip 
where the central gas pressure is 



m r , 



or /i 

P = P(0) = 2p A: B T /m p = — fiS 



r * 2 i 


4 








2 


L zl- 




^0 J 



(31) 




105 



128 y m r 




105 



p 

2 Mt /fo not 



STY 




Q!_ kB_ T\at_ 

fl V m p a S B - 

1/7 



1/7 



m„ cr S B 



(32) 



It is interesting to note a very weak dependence of Po on Q' and Bq, as well as on the 
opacity k, and a relatively strong dependence on E and fi. 



Likewise, the local plasma-/? parameter associated with the vertical field B is 

z 2 



P{z) = 8ttP/B 2 = ft 



1 - 



(33) 



with ft = 0(0) > 1 



/So = /3(0) 



P, 



\ 87T/ \/ m. 



P2/87T 128 

being the midplane value; our model assumes (3q 1. 




105 6 



sat 



P 



m p cr SB 



(34) 



This completes the calculation of the interior vertical structure of an optically thick 
MRI-heated gas-pressure-dominated accretion disk threaded by a relatively weak vertical 
magnetic field. 
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4. Breakdown of the Model near the Disk Photosphere 

The above model of the disk interior vertical structure relies on several assumptions that 
should be well justified deep within the disk but are expected to break down close to the 
disk surface, as z — > Zq. This situation is similar to the one in stellar structure calculations 
as one approaches the star photosphere. These assumptions can be cast in terms of a certain 
ordering of the relevant physical length scales which can be summarized as follows: 

Imn(z), X pb (z) < H{z) <Az<z . (35) 

Here, H(z) = —{dhiP/dz)~ l is the disk local gas pressure scale height, and Az = zo — z is 
the geometrical depth, i.e., the distance from the disk edge. 

In this section we shall estimate how rapidly each of these scales varies as one approaches 
the disk photosphere and will thus estimate at what optical depth the individual components 
of the above ordering break down. 

Here we are interested in the region near the disk's surface: 

Az = z — z <C z , (36) 

or, equivalently, 

AC = Az/z < 1 . (37) 

In what follows, it will often be convenient to use the optical depth coordinate r(z) = 
f*° np(z)dz, along with z itself; the two quantities are related to each other by equation (I2ip . 
In particular, near the disk edge, A( C 1, we have 

35 

t(z ~ zo) ~ 2k P (0)z AC 4 = — r tot A( 4 . (38) 

o 

Since equation ( )2~Tj) was derived for the optically thick region, the above expression is valid 
only for r > 1, i.e., for 1 > A( > 0.7r tot 1/4 . 

It is also interesting to note that because the gravity is nearly constant near the disk 
edge, g z (z ~ z ) = —Q 2 z ~ — Q 2 z = const, the vertical hydrostatic pressure balance, 
equation, eq. (JSJ), yields 

P(t) -r\g z (z )\/K^rn 2 z /K. (39) 

That is, the variation of the optical depth r with the geometrical depth Az is the same as 
that of the gas pressure in this region, and hence r can be used as a proxy for the pressure 



c.f. (e.g.. IShakura fc Sunyaevlll973t iHirose et al.ll2006l ). 



- 19 - 



Now we are in a position to check the validity of some of the assumptions underlying 
our model. 

(i) Gas scale-height, H. First, for our model to be valid near the disk surface, we 
must require that Az = z — z > H(z). 



The gas pressure scale- height, using equation (13*11) . is 

Zq 1_ 

c 



Near disk edge, A( <C 1, this becomes 

H(z)~~Az, (41) 

so that the assumption H(z) < Az is marginally satisfied. 

(ii) Photon mean-free path, A p h- Similarly, because of the radiative diffusion ap- 
proximation, our model is valid only as long as the photon mean free path is sufficiently 
short, i.e., A p h <C H, Az. Using the above expressions, we have 

A P h = - = (1 - C 2 r 3 - ^ AC" 3 , (42) 

where A p h,o is the photon mean-free path at the disk midplane: 

1 16 

Aph ' = M0^ = 35"° Tt ° tl - (43) 
From this we see that near the disk edge 

V_fV_Vo a<t 4 -— r- 1 Ar 4 -— (u) 

H ~ Az ~ 2z AC " 35 tot C " r(z) ' 1 J 

Thus we see that that the condition that A p h(z) <C H(z) < Az is automatically satisfied as 
long as we are in the optically thick part of the disk, r(z) S> 1. 

(iii) MRI scale, / mri . Next, we want to check that Z mri is less than H and Az. We 
have: 

V A ( p\ ~ 1/2 ( 2 \ " 3/2 

^mri = ~TT = ^mri,0 ( ) = ^mri,0 ( 1 ~ C ) j (45) 



where 



^vPo 



is the MRI scale at z — 0. 
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Near the disk edge we then have 

/ mr . ~ AC" 3/2 = AC 3/2 , (47) 

or 

Ui * ^ (I f) ~ 0.31 z ^ f~ 3 / 8 , (48) 



2x/2 V35 

where f(z) = r(z)/r tot . 
Thus, we have 



~ (2/3Q)" 1 / 2 AC~ 5/2 ~ 1.8 (3~ 1/2 f^ 8 , (49) 



H(z) Az 



and so the condition l mT i(z) <C #(2) is satisfied as long as A( 3> (2/3 ) 1//5 or, equivalently, 
f >2.5/3 " 4/5 . 

(iv) Finally, the calculation presented in the preceding section neglects the magnetic 
contribution in the vertical pressure balance, which is valid only as long as the local plasma-/? 
is greater than 1. According to equation ( I33p . in the outer layer of the disk, the plasma- 
parameter can be written as 

P(z -+ z ) = A, (1 - C 2 ) 4 * 16A) AC 4 ^ ^ /So f . (50) 

35 

Thus, the assumption (3 3> 1 is valid only for f 3> (35/128) /Jq" 1 . At optical depths below 
this, our model for the disk structure becomes invalid. Notice, however, that for /3 ^> 1, 
the critical optical depth corresponding to j3 = 1, f = (128/35) /Jg" 1 , is smaller than the 
critical optical depth f = 2.5 (3q^ 5 at which Z mri becomes equal to the gas scale height H. 
This implies the condition Z mr ; < H is more restrictive, and hence more important, than the 
condition j3 >> 1. 



5. Parker Instability and Coronal Power 



Let us now consider the Parker instability in the disk, with an ultimate goal of estimating 
the coronal fraction of the accretion power. We wish to remind the reader that the spirit of 
our approach is to consider the Parker instability as a secondary parasitic instability fe eding 
on the horizontal magnetic field of the primary MRI mode (e.g., iTout fc Pringldll992l ). and 
competing with the usual GX94 parasitic instabilities. Our approach is to compare the linear 
growth-rates of these two types of parasitic modes and see under what conditions the Parker 
instability becomes important. 
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The linear growth rate of a maximally-unstable parasitic instability scales as 7gx ~ fife, 
where b = B^ OT / B is the MRI channel flow 's horizontal magnetic field normalized by the 



vertical magnetic field ( IGoodman fc Xulll994l ). The growth rate of the fastest growing Parker 
instability can be estimated as 

1p{z) ~ V Ahoi /H(z) = bV A /H = Qbl mTi /H. (51) 

Thus, the ratio of the growth rates of the two types of parasitic instabilties is just proportional 
to the ratio of the MRI scale Z mri to the disk local pressure scale height H(z): 

1p k 



lax H 



(52) 



As we see, this ratio is small in most of the disk's volume and hence the GX94 parasitic 
instabilities win over Parker. However, as one approaches the disk surface the above ratio 
increases because H (z) starts to decrease due to decreasing temperature and increasing 
gravity, while, simultaneously, the MRI scale l mri ~ B Q /fly/A7rp increases due decreasing 
density. In particular, using our estimate (T4"9"j) near the disk surface, we get 

TP ^mri n-1/2 5/8 

IT T 1 . (5dJ 

1GX ti 

Thus, at a characteristic optical depth of order 

r P = r tot (3- 4/b (54) 

the local MRI scale / mr j becomes comparable with H and hence the the growth rate of the 
secondary Parker instability becomes comparable with that of the GX94 parasitic instabil- 
ities. We conjecture that, in the region above the critical depth zp = z(rp), the Parker 
instability starts to compete with other parasitic modes as the main mechanism of destroy- 
ing MRI channel flows. At the same time, however, we acknowledge that, since the condition 
A mri <C H is no longer satisfied in this region, the cla ssical incompressib le analysis of GX94 is 



no longer valid, and the effects of stratification (e.g.. lLatter et al.ll2010l ) need to be taken into 
account. Nevertheless, it is important to note that since, according to ( |50|) . the plasma- j3 
parameter at r = t p is still greater than 1 (it is of order (3 P ~ f3^ 5 3> 1), there is still a 
sizable region above zp where MRI, although modified by stratification, continues to operate 
and leads to release of the gravitational energy. It then seems reasonable to suggest that a 
substantial fraction of the accretion energy released at heights above zp is transported buoy- 
antly into the corona. Thus, a reasonable estimate for the thickness Azp of this buoyantly 
active region is given by the corresponding local MRI length-scale (which is of the order of 
the local pressure scale height at r = r P ), which can be estimated using eqs. (1451) and (154"|) : 

Az P ~ Zmri (tp) ~ z /3q 1/5 < z . (55) 
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Since the reduced volumetric MRI energy-release rate Q' is roughly uniform across the 
disk (see § !3.ip . we come to the conclusion that the fraction / of the accretion power released 
in the corona should scale simply as 

/ ~ Az P /z ~ P~ 1/5 . (56) 

Using equation ( |34"|) . the dependence of the coronal fraction on the primary input pa- 
rameters of our model can be expressed as 

/ocQ-^E-^(g) 6/35 r if . (57) 

Thus, within the limitations of our theory, the coronal power fraction is rather insensitive 
to most of the input parameters. 

Recall now that one of the assumptions of our model was that the vertical net field Bq, 
while small, is still larger than the field that would be generated by the large-scale MHD 
dynamo in the zero-net flux case. This limits the midplane plasma-/? parameter to something 

— 1/5 

of order 100, with the corresponding coronal fraction no less than / min ~ f3 Q max ~ 0.4. 



6. Mass Accretion Rate 



According to IShakura fc Sunyaevl (119731 ) , the mass accretion rate M is related to the 
total energy dissipation in a Keplerian disk via 



2 / Q'dz = 2 Q'z = — MQ 2 

J 87T 





(58) 



This allows us to express the mass accretion rate in terms of the local quantities in our 
model: 



M 



64tt Blhl 



sat 



3 87rfi 2 



c s ,o ■ 



(59) 



Substituting our expression ( )28|) for c Sj o we get a relationship between M and £>o, ^, 
and S (or r tot ): 

, r>2r9 8/7 /i*„\4/7 /7 - N 1/7 

M[B (R)MR)MR)]*tt(v^) fi~ 2 (— ) (— ) • (60) 

V 07T / \Tflp/ ^C^SB' 

As we see, the mass accretion rate has a relatively strong dependence on the magnetic field 
and on the rotation rate, but a rather weak dependence on the total disk column density 
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and optical depth. This effective insensitivity of M to r tot can be exploited. 

1 /7 

ignoring the r tot dependence and taking the profile to be Keplerian, fix ^ 



For example, 
R~ 3/2 , we get 

a scaling M ~ (R) R 3 . Thus, we see that a stationary accretion regime, M(R) = const 
requi res a special particular mag netic field profile, B-Q t£lt (R) ~ i? -21 / 16 , which is in agreement 
with IShakura fc Sunyaevl (119731 ) . 



In general, equation (j60l) gives us an equation that governs a longer-term evolution for 
the gas surface mass density: S(-R, t) = dM /27iRdR, relating the evolution of the disk density 
to the radial distribution of the vertical magnetic flux through the disk. Determining the lat- 



from being solved (e.g. 


van Balleeooiien 


Il989: Lubow et al. 1994: Soruit & Uzdenskv 2005: 


Uzdenskv & Goodmanl2008: 


Rothstein & Lovelace 2008; Lovelace et al.ll2009l; Beckwith et al. 


2009; Fromane & Stone 


200< 


))• 



Vertical Structure of the Disk with Zero Net Flux 



The picture presented in the previous sections was developed under the assumption 
that the super-imposed net vertical magnetic field, Bq, while weak compared with the gas 
pressure inside the disk, nevertheless dominates over the field -Bdyn that would be generated 
by the MRI-driven turbulent dynamo in the absence of B . Since, as numerical simulations 
demonstrate, Bd yn is t ypically indeed rela t ively weak, corresponding to a cen tral plasma-/? of 
order /3d yn ~ 100 (e.g., iHirose et al.ll2006t iDavis et al.ll2010t IShi et al.l 120101 ). there is indeed 
a sizable range of parameters where the condition B\ jn <C Bq -C B^ q is satisfied and hence 
where the above picture applies. 

However, it is also interesting to consider the zero net flux case where there is no ex- 
ternally imposed vertical magnetic field (or where this field is weak compared with -Bdyn)- 
This case has attracted a lot of attention in the MRI literature, especially in recent years (e.g., 



Hirose et al 



2010 



2006 



Fromang &: Papaloizou 



Guan fc Gammie 



2011 



2007: 



essah et al. 2007: Davis et all2010: Shi et al 



Simon et al.ll201l[ ). In this case the magnetic field responsi- 



ble for driving the MRI is the self-generated turbulent field produced by an MHD dynamo 
associated with the turbulence itself. One expects both small-scale and large-scale dynamo 
generating the field on a broad range of spatial scales. The problem of computing analyt- 
ically the overall spectrum of the resulting magnetic field and its vertical distribution in 
a stratified disk is an outstanding problem in accretion disk theory. This formidable chal- 
lenge has so far eluded a f ull solution, although several important theoretical advance s have 



already been m ade (e.g., iLesur fc Ogilvid l2008t IVishniad l2009t iBlackman fc Pessahl 12009 



Blackmanl 120121 ) . It is, however, reasonable to expect that the characteristic length scale of 
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the average magnetic pressure (averaged over turbulent fluctuations) is much larger than the 
dominant MRI wavelength Z mri and, especially if a strong large-scale MRI dynamo is active, 
may even be larger than the gas pressure scale height H. In that case, in line with the ar- 
guments presented in § 13.11 we may expect the volumetric dissipation rate Q(z) ~ QB 2 /8ir 
to be also roughly uniform in z within the disk. We note that this conclusion differs dras- 
tic ally from the assumption t hat Q(z) traces the mass density profile, Q ~ p(z), adopted 
by lShakura &: Sunyaevl (119731 ). However, our c onclu sion is in a good agre ement with the re- 
sults of numerical simulations by lHirose et al.l ( 120061 ) and lShi et al.l (120101 ) — the only known 
to us numerical studies of MRI turbulence in a gas-pressure-dominated disk with optically- 
thick radiative cooling — who report a flat top-hat energy dissipation profile over 3 scale 
heights on each side of the dis k. [The flat magnet ic energy and dissipation profiles are also 
consistent with the findings by ISimon et al.l (120111 ) for isothermal stratified disks.] If this is 
indeed the case, then we can apply the rest of our analysis presented in § 13.21 to the zero- 
net flux case, except that Bq/8tt should be replaced with B^ yn /87r ~ P . In particular, 
we still recover the z-profiles of p and T as those given by equations (jT9l) and (|2~TT) . e.g., 
a parabolic profile for the temperature T(z) = T (1 — z 2 /zq). We find that these profiles 
agree well with the actual average profiles measured in numerical simulations by lHirose et al. 
J2006I ) (see their Figs. 2 and 3). 



It is important to note that if the zero-net-flux case the characteristic magnetic field 
strength can no longer be viewed as an external input parameter but rather has to be 
determined self-consistently along with the other disk quantities. The only input parameters 
in this case are Q and E, plus the opacity k (which we assume to be constant, as is the 
case of the dominant electron scattering). In lieu of using the magnetic field as an input 
parameter, we can just use the result /3q = (3d ya ~ 100 obtained from numerical simulations. 
Substituting this on the left-hand side of equation (|34|) for (3 , we can rearrange it to express 
the characteristic magnetic energy density in terms of Q, E, and k, and then substitute the 
resulting expression for the magnetic energy into our expressions (12"B1) - (13"01) for T , z , p , etc. 
Ignoring factors of order unity we get: 
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We can then also use equation ( 1601) to express the mass accretion rate in terms of £ 
and Q: 



M 



M 4/3 (^VV^ir^/s 
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(66) 
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As one can ea s ily ch eck these scalings are exactly the same as those obtained by 
Shakura &; Sunyaevl (119731 ) [see their eqs. (2.16)] in the appropriate limit (gas-pressure dom- 
inated disk, constant opacity, etc.). 

Thus, we believe that much of the analysis developed in the main part of this paper 
can also be applied to the zero-net-flux case, resulting in an essentially very similar vertical 
disk structure (repr esented by the temperatur e and density profiles) and enabling one to 
recover the classical IShakura fc Sunyaevl ( 19731 ) scalings for the key disk parameters. While 
in this paper we do not want to make any strong claims regarding the coronal fraction / of 
the accretion power in the zero-net-flux case, we nevertheless note that a typical value of 
0o ~ 100 expected in this case implies, in conjunction with equation ( 156|) . a relatively high 
(tens of percent, which is consistent with AGN observations) and universal (independent of 
any system parameters) value of /. 



8. Conclusions 

In this paper we considered the problem of the vertical thermal structure of a thin gas- 
pressure-dominated accretion disk heated by the dissipation of MRI turbulence and cooled 
by optically thick radiative cooling. We also considered the question of the fraction of the 
overall accretion energy that is transported by buoyantly rising magnetic loops into the 
tenuous corona lying above the disk. 

In the main part of the paper, we focused on the non-zero-net vertical flux case, in 
which the disk is threaded by an externally imposed vertical magnetic field, B . Because 
we are interested in an MRI-active disk, we considered the case where this field is relatively 
weak, Bq/8tt <C Pq, where Pq is the gas pressure in the middle of the disk. At the same 
time, we assumed that the mean vertical field Bq is stronger than the dynamo-generated 
magnetic field Bd yn , which, according to numerical simulations, is expected to be of order 
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B^ yQ /8iT ~ 10 2 Pq. We argued that, under these assumptions, the volumetric heating rate, 
Q, due to the dissipation of the MRI turbulence (which w e tentatively associate with the 



disruption of MRI channel modes by parasitic instabilities, iGoodman fc Xul (119941 )) should 
scale as QB 2 /8ii] in particular, it should be independent of the height z inside most of the 
disk. It is particularly important that the ^-profile of Q does not trace the profiles of the 
gas density or pressure. Making use of this finding and assuming, in addition, that the 
opacity k is also constant in z, we then were able to solve analytically the combined set of 
equations governing the vertical structure of the disk — the hydrostatic pressure balance, 
the energy conservation, and the optically thick radiative transfer equation. As a result, we 
were able to obtain the z-profiles of the the gas temperature and density inside the disk: 
T(z) = T (1 — z 2 /zq), p(z) = po (1 — z 2 / z o) 3 , where T and po are the values at the disk 
midplane, z — 0, and z is the effective thickness of the disk [see eqs. ffTUj) and (j2"T|) ]. 

We were also able to to determine all the key disk parameters, such as the midplane 
temperature, density, and pressure, and the disk thickness zq, in terms of the governing input 
parameters in this problem: the external vertical magnetic field Bo, the surface density S, 
the disk rotation rate Q, and the opacity k. This enabled us to evaluate the scaling of the 
mass accretion rate M with these parameters (see § E]) and formulate the equation governing 
the time evolution of the radial distribution of mass across the disk, £(r, t). 

We then also examined how and where various assumptions on which our disk model is 
based — (3(z) 3> 1, / p h <C H , etc. — break down as one approaches the disk's surface. This 
enabled us to address the question of the coronal power fraction. In our view, the coronal 
power is governed by the competition between various parasitic instabilities disrupting the 
primary MRI channel modes: the GX94 instabilities leading to fully developed local MHD 
turbulence whose dissipation heats the disk locally, and the Parker instability that pumps 
the Poynting flux into the ADC by buoyantly rising magnetic loops. The magnetic energy 
of these loops is then dissipated by reconnection in the corona. In our model, the resulting 
corona fraction turns out to be relatively insensitive to most input parameters, scaling as 
/ ~ A) 1 ^ ■ The practical consequence of this conclusion is that it is difficult to avoid a sizable 
(tens of percent) fraction of the accretion power to be released in the overlying corona. The 
part of this dissipated energy that goes to the electrons powers the observed coronal emission, 
whereas a significant part of the ion coronal energy may actually be transported back to the 
disk by ions streaming along the closed field lines and get deposited in the dense disk through 
ion-ion collisions. Finally, some of the energy dissipated by reconnection involving open field 
lines may power outflows (winds and jets) along these open field lines. 

Finally, in § 0, we turned our attention to the case of a disk with a zero net flux. This 
is the case where the next externally imposed vertical field Bq is either absent altogether or 



-27- 



small compared to the MRI dynamo-generated field -Bdyn- In this case, the entire steady-state 
vertical structure of a thin gravitationally-stratified disk should be determined solely by S, 
Q, and k since there are no other parameters in the problem (ignoring external irradiation). 
As numerical simulations show, in this case one expects an effective large-scale MRI dynamo 
to produce a sizable large-scale (comparable or larger than the disk pressure scale-height H) 
magnetic field B\ yn /&-K ~ f3^ n Pq with /3d yn — 10 2 . We argued that since the magnetic energy 
density is roughly uniform in z inside the disk (perhaps up to the equipartition height at 
which P(z) = B^ yn /8ir ~ 10~ 2 Po)> then the volumetric energy dissipation rate of the MHD 
turbulence should also be roughly uniform, just as it is for the finite net flux case. This 



pictu r e is in fact support ed by numerical simulations of stratified disk (e.g.. iMiller fc Stone 



20001 ; iHirose et al.ll2006l ). We then argued that we can apply the analysis developed in 
§ 13.21 to the zero-net-flux case. As a result, we recover the same vertical profiles for the 
gas temperature and density, i.e., T(z) = T (1 — z 2 /zq), and p(z) = p (1 — z 2 /zq) 3 [see 
equations (fT9l) and (l2~Tj) ]. The se profiles seem to b e a good agreement with those obtained in 



full numerical simulations by lHirose et al.l ( 120061 ). Furthermore, we find our scalings of the 



main accretion disk parameters (the midplane values of the temperature and density, the 
disk thic kness Zg, and the mas s accre tion rate M, etc.) to be the same as those in the classic 



paper by IShakura fc Sunyaevl (119731 ) for the regime under consideration (constant opacity, 



gas pressure-dominated disk). 

While the theoretical model presented in this paper is successful in being able to pro- 
vide a set of concrete, physically motivated predictions, it still has to be viewed just as an 
idealized conceptual toy model. It relies on a number of simplifying assumptions — e.g., gas- 
pressure domination, a single constant opacity, external irradiation — that preclude its direct 
application to real astrophysical systems. For example, the inner regions of optically-thick 
BH accretion disks are dominated by the radiation pressure, which is completely ignored in 
the present study. In addition, the magnetic field generated by the MRI dynamo is either 
completely ignored (in the first half of the paper, where we consider a disk threaded by an 
external magnetic field) or is treated in a simplified fashion (in §[7]). Perhaps for these reasons 
the present model is not able to address the questions of spectral state transitions in galactic 
BH binaries and, in particular, to explain the relatively low level of X-ray coronal activity in 
the high-soft state. Future theoretical studies should develop generalizations of the present 
work to take into account radiation pressure, more realistic opacities, and, perhaps, external 
irradiation. One can also envision a more rigorous and, perhaps, more accurate model of 
MRI saturation and MHD turbulent dissipation; the development of such a model should 
benefit from detailed direct comparisons with numerical simulations. 



In conclusion, we believe that the future of accretion disk studies lies in incorporation 
of more realistic radiation and thermal physics — i.e., more accurate treatment of the disk 
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thermodynamics and radiative cooling processes. We hope that the present paper will help 
stimulate and pave the wa y for such studies using a dvanced n u meric a l simulations , alon g 



the lines of the studies by iHirose et all ( 120061 . 120091 ): IShi et all ( 12010( 1: iBlaes et all ( 1201 lh 



Eventually, we hope, this line of research will reach the state of maturity at which it is able to 
provide meaningful predictions for observations and explain important observational facts. 

I am very grateful to Prof. Jeremy Goodman for many fruitful and stimulating dis- 
cussions. I also would like to thank M. Begelman, J. Goodman, and J. Simon for useful 
comments on this manuscript. This work is supported by National Science Foundation 
grant No. PHY-0821899 (PFC: Center for Magnetic Self-Organization in Laboratory and 
Astrophysical Plasmas) and by NASA under grant No. NNX11AE12G. 
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